Strong disorder effects of a Dirac fermion with a random vector field 
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We study a Dirac fermion model with a random vector field, especially paying attention to a 
strong disorder regime. Applying the bosonization techniques, we first derive an equivalent sine- 
Gordon model, and next average over the random vector field using the replica trick. The operator 
product expansion based on the replica action leads to scaling equations of the coupling constants 
( "fugacities" ) with nonlinear terms, if we take into account the fusion of the vertex operators. These 
equations are converted into a nonlinear diffusion equation known as the KPP equation. By the 
use of the asymptotic solution of the equation, we calculate the density of state, the generalized 
inverse participation ratios, and their spatial correlations. We show that results known so far are all 
derived in a unified way from the point of view of the renormalization group. Remarkably, it turns 
out that the scaling exponent obtained in this paper reproduces the recent numerical calculations 
of the density correlation function. This implies that the freezing transition has actually revealed 
itself in such calculations. 
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I. INTRODUCTION 

A Dirac fermion model has proved to be quite useful, 
even in condensed matter physics, as an effective theory 
for various kinds ofjphcnomena, e.g, the integeE-o 
turn Hall transition,Em'[j cZ-wave superconductorsJJt 
etc. It is also involved with statistical mechanical models 
such as the Ising models or the XY-modelJlJ Especially, 
including randomness, this model has clarified remark- 
able aspects of disordered systems, for example, lnlnr 
behavior of the specific heat of the randompb«&dJsing 
modellli and multifractal scaling dimensions .ErH't3EJ 

Many problems, however, still remain to be explored. 
Especially, the critical theory of a Dirac fermion with 
generic disorder, which is believed to be an effective 
theory of the integer quantum Hall transition, is still 
missing.B In this case, coupling constants describing dis- 
order strength flow to a strong coupling regime, which we 
cannot reach by weak coupling approaches. Even simpler 
model with a random vector field only, whose zero energy 
state can be obtained exactly for any realization of ran- 
domness, yields non-trivial strong disorder phase, recog- 
nized through the calculatian-of the generalized inverse 
participation ratios (IPR) OB 

On the other hand, Carpentier and Le Doussal have 
recently developed a—theory for the XY-model with a 
random gauge field. Ej It has been shown that a naive 
replica theory, for this model concludes a reentrant 
transitionJiaO which is not consistent to the numerical 
calculationseZI pad to exact arguments by the use of the 
Nishimori-line.Efl Although this inconsistency has been a 
long standing jjamjplem, underlying physics has recently 
been clarified o£J In particular, ScheidlE] has pointed 
out the importance of a role which vortices with higher 
charges play. Carpentier and Le Doussal have derived 
new scaling equations for the Coulomb gas model tak- 
ing into account "fusion" of vortices, from which higher 



charge vortices are generated consistently. The renor- 
malization group (RG) equations they have derived have 
intimate relationship with the freezing transition of spin 
glasses (the random energy model) .Ell 

By the use of the fact that the exact zero energy wave 
function of a Dirac fermion model with a random vector 
field is equivalent to the vertex operator of the boson field 
theory,E3 Carpentier and Le Doussal have applied their 
RG method to a Dirac fermion model, calculating the 
IPR. Variational method has also been .applied in order 
to calculate the density of state (DOS). £21 It has turned 
out that these quantities shows peculiar behavior in a 
strong coupling regime due to the freezing transition. 

Motivated by these developments, we rederive, in this 
paper, scaling equations for a replicated sine-Gordon 
model with a random vector field, which is the bosonizcd 
model of a Dirac fermion with a random vector field. The 
main point is that we take into account the fusion of the 
vertex operators. In order to derive scaling equations 
at one-loop order, we use the operator product expan- 
sion (OPE) techniques. Based on these, we calculate the 
DOS, the IPR, and their spatial correlations of a Dirac 
fermion model. We show that it is possible to derive all 
results known so far in a unified way from the RG point 
of view. Moreover, it turns out that the scaling exponent 
obtained in this paper reproduces the recent numerical 
calculations tw-|Ryu and Hatsugai of the density correla- 
tion function.EJ This implies that the freezing transition 
has manifestly revealed itself in their calculations. 

This paper is organized as follows: In the next section, 
a Dirac fermion model is introduced and a correspond- 
ing sine-Gordon model is derived. In Sec. Ill , the OPEs 
are calculated, which lead to scaling equations, by taking 
account of the fusion of the vertex operators, and in Sec. 
IV , they are converted into more convenient form, known 
as the KPP equation.OE3c3 By the use of the asymptotic 
solution of the equation, the DOS and the IPR are cal- 
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culated in Sec. A summary and discussions are given 
in Sec. Ivi. 



II. DIRAC FERMION WITH RANDOM VECTOR 
FIELD 



In this section, we first introduce a Dirac fermion 
model including a random vector field and give formulae 
to calculate the DOS, etc. Next, using the bosonization 
techniques, we convert the Dirac action into an equivalent 
sine-Gordon action. The latter is the model we directly 
investigate in this paper. 



A. Dirac fermion model 

The action functional is defined by 

S= j d 2 x (d„ - iAJ i/j - ijjW] , (2.1) 

where y = u — iE, 7 M = ay, (the Pauli matrices), and 
denotes a quenched random vector field with the follow- 
ing gaussian probability distribution 



P[Afj] (x exp 



1 

2tt5 



d 2 xA 2 



(2.2) 



Here we include an extra factor 7r into the definition of 
the distribution width g for later convenience. The Green 
function is computed as 



Im tv(ip(x)ip(x)} 
= tr(x 



[z 7M (cV-a M )-£:] 2 +^ 



(2.3) 



where tr is the trace with respect to the two-component 
spinor. Let us further define the qth power of the Green 
function 

T (q) (x) = (Im tr(^(s)^(ic)»* . (2.4) 
Then, the following relation is valid 
uj q ~ l T {q \x) 

^C q J2\*n(x)\ 2q S(E n -E), (w ->+()), (2.5) 



where ^ n {x) is the eigenstate of the Hamiltonian with 
energy E = E n , namely, ij^d^ - iA^n = E n H? n , and 
C q is a numerical constant C q — Tr(2q — 3)!!/(2q — 2)!!. 
Field-theoretically, the ensemble-averaged DOS p(E), 
the generalized IPR P^(E), and thpiii-.spatial corre- 
lations Q {quq2 \x,y,E) are defined by&Ea, in the limit 



p(E) = rm( x )/c 1 



P^(E)=Lo q - 1 r(i)(x)/{C q p(E)}, 
Q^^\x~y,E) 

= u J q ^- 1 T^)(x)TM(y)/ [C qi C q2 p{E)\ . (2.6) 

Here we have assumed that the ensemble-average recovers 
the translational invariance and hence T^(x) does not 
depend on x and r^i) (x)T^ 2 '> (y) depends on the separa- 
tion x — y only. 

B. Sine-Gordon model 



Bosonization of (^lj) yields the following sine-Gordon 
action 



fx 

47T 



{d^cj)) 2 + 2iA ll e l _ lv d v (j) - y cos <j> 



(2.7) 



By the use of the replica trick, ensemble-average over the 
vector field A^ gives 



Z m 



V(f>e 



(2.8) 



with 



ax 

47T 



- ^2 d^gGgbd^b - y cos <j> a 



a,b=l 



o=l 



where a, b denote m replicas and 



G a b 



1 

K 



Sab + 9- 



(2.9) 



(2.10) 



The bare coupling constant is given by K = 1, since the 
present system is a free fermion model. 

So far we have derived a replicated sine-Gordon model 
~9|) by the bosonization of the Dirac fermion model. 
We also reach the same model from the XY-model with 
a random gauge field as a result of applying the Villain 
approximation, "integrating" over vortices, and_fj>erform- 
ing the quenched average with a replica trick. tj In this 
case, the coupling constant K is K = J/T, where J and 
T denotes, respectively, the spin coupling constant and 
the temperature. The problem had been that the RG 
analysis based on this action (2.£) predicts apceentrant 
transition, which contradicts numerical studiestland ex- 
act arguments by the usf of .the Nishimori-line.lla 

Recent development sJiSE3E-a however, have made it 
possible to describe the random XY model within a 
framework of the replica method, without resort to the 
replica symmetry breaking. In the language of the sine- 
Gordon model, the main point is that higher charge ver- 
tex operators, which are not included in the above action 
but are generated in the process of the renormalization, 
are taken into account. The next section is devoted to 
the reformulation of the method. 
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III. SCALING EQUATIONS 

Carpentier and Le Doussal have recently proposed iro-, 
proved scaling equations for the Coulomb gas model.tll 
They have taken into account the fusion of vortices, 
which results in scaling equations with nonlinear terms. 
Without these terms, the flow of the coupling constant y 
(also called "fugacity" in the context of the Coulomb gas 
model) is governed by the expo nent 2 — x, where x is the 
dimension of ^ a cos <j) a in Eq. (2.9). The nonlinearity of 
the scaling equations, however, has turned out to play a 
crucial role in the RG flow. In what follows, we rederive 
the scaling equations for the present sine-Gordon model 
by the use of the OPE techniques. 

We start with the same sine-Gordon action but with 
generalized fugacities, 



5* 



fx 

ill 



-d^aGabd^b - 22 Y(n)e 



(3.1) 



where a, b — 1,2, •••to, n* = (m, 112, ■ • ■ , n m ) is a 
to component vector with integer elements n a which 
will be specified momentarily. The vector n denotes a 
"charge" of the vertex operator, so that let us define a 
total charge n = ^2 a n a . Initially, the bare fugacity is 
Y(n) = 2/i(5 n ,m + <5 n ,-ni), where the vector ni is such 
that only one of its elements is 1 and others are zero. In 
this case, X) n #o Y{n)e m ^ a = 2y x Y%=\ cos (0a), and the 
action fl3.1|) reduces to the normal replicated sine-Gordon 
action (|2.9[). 

The unperturbed correlation function reads 



= G- b 1 ln\x-y\- 2 . 
The inverse of the matrix G is easy to calculate, 

9K 2 



(3.2) 



G ab = KS ab - " 



mgK 



KS„ 



9K 2 



(3.3) 
in the last 



where we have taken the replica limit to 
equation. By using the correlation function (3.2), we 
readily find, 



i<t>a(z) e -in a 4> a (0)\ _ 



where the dimension x(n) of e ma ^ a 
x(n)= n'G" 1 !! 



(3.4) 



is 



= K\n\ 2 - gK 2 n 2 , 



(3.5) 



with |n| 2 = n • n and n • n' = J2 a n a n ' a - 

The OPE of the vertex operators in the case n + n' = 
n" ^ is given by 



e in a <t>a(z) e in' a <j> a (0) _ 



|^|x(n,n') 



in"<p a (0) 



(3.6) 



with 



x(n, n')= x(n) + x(n') — x(n") 
= -2n t G" 1 n' 
= 2 (gK 2 nri ~ Kn ■ n') , 



(3.7) 



where we have used the fact that the matrix G is sym- 
metric. On the other hand, in the case n + n' = 0, we 
find 



1 



e 'e 



|^|2x(n) 

-1 
\ z \2x(n)-2 ( 



6 a (nn*) afc a^, (3.8) 



where z, z = x± ± ix2, d — d z and d = d s 



First, the OPE (|3.6| ) leads us directly to the scaling 
equations, 



(3.9) 



where I = In L with the system size L. 

So far we have not specified the vector charges n. The 
initial condition at / = is such that there exist ni-type 
vector s on ly. However, the "fusion" of the vertex oper- 
ators (3.6) yields higher charge vectors in the process of 
the renormalization. Let us fix the charge n, and com- 
pute the scaling dimensions of the vertex operators. Con- 
sider the case with n = 1 for simplicity. Then, we have 
n* = (1,0,0, •••), (1,1, -1,0 •••), and (2, -1, 0, • • •), for 
example, whose dimensions are, respectively, K — gK 2 , 
3K — gK 2 , and hK — gK 2 . It is readily seen that the vec- 
tors of the first type have the most relevant dimension. 
In a similar way, given a charge n, we find that the most 
relevant vectors are the ones that have |n| l's (— l's) for 
a positive (negative) n. 

In what follows, we restrict ourselves to these most 
relevant ones, n with n a = 0, 1 and with n a = 0,-1, 
and therefore, vectors whose charges lie between ±m are 
taken into account. Let us define a set of vectors N g 
which includes vectors n with n a — 0, 1 and ^ a n a = q. 
The scaling dimensions depend, in this case, only on |n|, 
so that let us define, for n > 



x(n) = x n = Kn — gK 2 n 2 . 



(3.10) 



Then the scaling equations (3J)) decouple into positive 
and negative charge sectors, which obey the same scal- 
ing equations with the same initial condition. Thus, it 
turns out that Y(— n) = Y(n). Moreover, Y(n) depends 
on n only. Denoting Y(n) = Y{— n) = y n and counting 
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t he multiplicity in the summation over n' and n" in Eq. 
( |3.9|) , we end up with a closed set of scaling equations 



dl 



1 ™ f n\ 

(2 - Kn + gK 2 n 2 ) y n + - ^ ( ^jy^Vn-n'- 

(3.11) 



4 ' vn' 

n' = l 



Next consider the OPE (3.8). The exponent in this 
equation is 2x n — 2 = 2(n — gn 2 — 1), where we have 
used the free fermion condition K = 1. It should be 
noted that in a weak disorder regime, a naively repli- 
cated model including y\ only could describe the crit- 
ical behavior correctly. The exponent of this fusion is 
2x\ — 2 = —2g, which is definitely negative, and there- 
fore, gives rise to no ultraviolet singularity. Hence, the 
kinetic term should not be renormalized. This may be 
valid in a strong coupling regime, since large g tends to 
make the exponent negative, although in some range of n, 
positive exponents appear if one includes higher charge 
vectors. Taking these into account, we postulate that the 
kinetic term is not renormalized in any regime of interest 

K ~ 1, in what follows. 

The Green function can be described by the replicated 
bose fields as r''' ~ (cos</>i)(cos<?!>2) ■ • • (cos</> 9 ). If one 
keeps the most relevant field, one finds 



r(9>(a:) ~ (cOs(n a Mx))), 



r(9i)(x)rte)( y ) 

~ (cos (ri a <j> a (x)) cos (n't foiv))}, 



(3.12) 



after the ensemble average, where n £ N g , n' £ N gi , and 
n"eN 09 . 



IV. THE KPP EQUATION 

So far we have derived the scaling equations with 
nonlinear terms due to the fusion of the vertex opera- 
tors. Remarkably, it has been showntil that the non- 
linear terms are involved with the freezing transition 
of the random energy model. Ell This section is de voted 
to the review of converting the scaling equations (3.11) 
into the KPP e*Hiation,L§t3tII according to Carpentier 
and Le DoussalJi-J First, make the Mellin transformation 
2M(0 = Jo°° dzz n ~ 1 P(z, I), which is converted, by the use 
of the variable z = e* 
present purpose as 



into more convenient form for the 



y n (l) = / due nu P(u,l). 

J — oo 

In the w-space, n can be replaced by n 
we have 



diP(u, 1) =[2 + Kd u 



gK 2 d 2 u ] P(u,l)- 

O 

du^b (u — ln(e^ 



(4.1) 



-d v , and then 



-y P(u,l) 



xP(u u l)P(u 2 ,l). 



(4.2) 



Integration over u of the both sides of the equation yields 

1 



(4.3) 



It turns out that yo, which is just the normalization of 
the function P(u,l), flows as y — > 8. Therefore, setting 
yo = 8 gives 

diP(u, I) = [Kd u + gK 2 d 2 ] P{u, I) - 2P(u, I) 



— oo J —oo 



+- / d Ul I du 2 S(u- ln(e M1 +e" 2 )) 



xP(ui,l)P(u 2 ,l). 



(4.4) 



Although it may be difficult to solve this equation an- 
alytically, a crucial observation was made by Carpentier 
and Le Doussal as follows:li-3 Define a generating func- 
tional of y n by 



G(x,l)=l-(cxp(~ze- x+Kl )) p 



n+l 



(4.5) 



n=l 



where 



e u , and we have introduced the notation 



(A)p = ^P(u,l)A, since the function P(u,l) can 
be interpreted, after normalized to. unity, as the proba- 
bility distribution of the fugacity.E-1 What is remarkable 
is that the function G(x, I) obeys the following equation 



~d l G = D8£G + G(l-G), 



(4.6) 



where D ^.mKI/2. This equation is known as the KPP 
equation Fir tl3 Since yo = 8, the boundary condition is 
G(-oc, I) = 1 and G(oo, I) = 0. 



V. THE DOS AND THE IPR 



In order to calcu late the dominant exponents of the 
DOS, etc in Eq. (2.6), we need to know the scaling 
exponents of the local comp osite operators r( g ). Basi- 



cally they are given by ( 3.10 ) through the relation ( 3.12 ). 
However, it will turn out that the nonlinear terms of the 
scaling equations induce nontrivial dynamical scaling ex- 
ponents in a strong disorder regime. Therefore, we first 
have to read the dynamical scaling exponents from the 
solution of the KPP equation, and next calculate various 
exponents of the DOS, etc. 
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A. Dynamical scaling exponents of the vertex 
operators 

First of all, let us summarize the remarkable selec- 
tion rule of the front velocity of the KPP equation 
(4.6). Provided that G(x,l) has a boundary condition 
G(— oo, I) = 1 and G(oo, I) — as well as an initial con- 
dition G(x, I = 0) ~ e~^ x at large x. In this case, the 
KPP equation admits the following traveling solutionES 



G(x, 1) = g(x — m(l)), for I 



(5.1) 



The solution tells that the wavefront (a kink) travels with 
a speed dim(l) if we interpret / as time. Neglecting the 
nonlinear term G 2 and assuming the form of the solution 
as exp [— /i(i£ — vl)], we expect that the velocity v, defined 
by m(l) = vl, is given by v = 2(D/i + l//x). Actually, this 
is valid, as ikpg as y, is small enough. However, it has 
been shownc3 that if fi is larger than fi cr — 1/y/D, the 
relation does not hold any longer, and the velocity is a 
constant (the minimum _nalue) independent of fi. To be 
m(l) is given byn3 



precise, 



m(l) 



2{Dn + l/n)l-\ 
\fB_Ul- ±ln( 
£>(4Z- fin/' 



0(1) 
+ 0(1) 
+ 0(1) 



for 
for 
for 



n < i/VA 
n = i/Vd, 
n > i/Vd. 



(5.2) 



Next task is to obtain the behavior of the fugacities 
y n , w hich are re late d to the above solution of G(x, I) via 
Eqs. (4.1) and (4.5). To this end, notice that the front 
velocity .above is related to a typical value of u via the 
relation,oE3 

U t yp= [hiZ]p 

/OO 
dx [exp (-e- x+Kl ) - (exp (- ze~ x+Kl )) p ] 
-OO 

/OO 
dx [G(x, I) - {1 - exp (~e- x+Kl ) }] 
-OO 



m(l) - Kl. 



(5.3) 



Here, in the last line, we have used the fact that the wave- 
front of G(x, I) is located at x ~ m(l), whereas the wave 
front of 1 — exp (— e ~ x+Kl ^j at x ~ Kl. The distribution 
function P(u,l) is broad especially in the case of a pulled 
front (fi > 1/VD), and the typical value corresponds to 
the maximum of the distribution function. The typical 
Utyp may define typical values of y q ,typ through the rela- 
tion K%, 



Vq,typ 



\m(l)-Kl]q 



(5.4) 



Now we can define the dynamical scaling exponents z q 
from the /-dependence of the typical values y q ,typ{l)- To 
this end, let us recall the following fact. In a weak disor- 
der regime, the dynamical exponent of y\ is, as expected, 
just 2 — xi, where x\ is the scaling dimension of cosn a <?!> a 
with n 6 Ni. In a strong disorder regime, however, it has 



been shown£jE3 that the interaction terms in Eq. (3.1) 
with the initial condition y n (0) = yi$ n< i plays a crucial 
role in the RG flow, and therefore it is quite important 
to keep the operator y\ ^cosn a c/f> a in the action and to 
derive the scaling equation of y\ coupled together with 
those of higher y n generated by the fusion of the vertex 
operators. 

Similarly if one wants to know the dynamical scaling 
exponent of the higher y q , one has to perturb the ac- 
tion by adding the corresponding operator y q cosn a </> a 
(n £ N g ) as a source term. Namely, one has to consider 
the initial condition y n (0) — y q 5 n , q , and hence \i = q in 
Eq. fljT2| ), Thus, we find y q ,typ(l) ~ e 
defined by 



where z q is 



2„2 



Kq + qK^q 
l)Kq 



for 
for 



gK 2 < 2/q 2 , 
gK 2 > 2/q 2 . 



(5.5) 



Here we have neglected the logarithmic corrections ap- 
pearing in Eq. (5.2). The vertex operator cos(n a ^ Q ) with 



total charge q i.e, n € N g , therefore, obeys the scaling 
law with scaling exponent 2 — z q . 

In pa ssin g, we mention that the logarithmic corrections 
in Eq. (|5.2|) are universal and give corrections to the cor- 
relation functions like marginal perturbations, although 
we have neglected them in this section. We will briefly 
discuss them in the next Sec. [v|, calculating logarithmic 
corrections to the DOS. 



B. Calculation of the DOS and the IPR 

By using the dynamical scaling exponents, we first cal- 
culate the DOS defined in Eq. (£j]). Although we are 
basically interested in the Dirac fermion model, K = 1, 
we will derive formulas below using generic K, because 

they are in fact valid for K ~ 1. Notice that the energy 
E has the same dimension as y\. Then, we find 



A 

E 



e (2-K+gK*)l for gR 2 < 2) 

e (Vsg-iWi for gK 2 > 2i 



(5.6) 



where A is a renormalized energy, and in what follows, 
we set A = 1 for simplicity. On the other hand, since the 
scaling exponent of is 2 — zi, we have 



p{l) 



,-(K-gK 2 )l for gR 2 < 2j 

r (2 +f f-if^)l for gR 2 > 2 . 



This equation, together with Eq. (5.6), yields 

K-gK 2 

I E- 

P(E) 



E 2 - k +sk 2 for gK 2 < 2, 
E i-^-DK for g K 2 > 2. 



(5.7) 



(5.8) 



This result coincides with the one obtained by a varia- 
tional method.CJ In the weak coupling regime of the Dirac 
fermion, g < 2 with K = 1, this DOS is just the same as 
that obtained via usual weak coupling approaches B 
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Next, let us calculate the IPR according to the definition (|2.6D. The IPR are expected to show a power law behavior 
p(i) ~ L' T{q) . Since the energy E scales as ([[|), the IPR obey the scaling law P (q) {E) ~ E r(q)/z " as a function of 
the energy. One can read in this way the exponent r(q) from P^ q '{E) in what follows. In the case of gK 2 < 2, the 
qth moment scales as 



1 e -(q-l)(2-K+gK 2 )l-(2-^K q +Kq)l fo r g > ^2/gK 2 , 1 ' 



where we have used the fact that u> obeys the same scaling law as E in Eq. ( p.6|) . In the same way as the DOS, we 
reach 



{q-l){2-gK* q) 



P^E) ~ E Z^, " K < *)■ (5.10) 

[ £ 2-k+ 9 k^ for g > ^/2/gK 2 . 

Contrary to this, in the case gK 2 > 2, we have 

u*-im fe-^t^-^H^V)! fOT q< ^2/gjO, 



e 



-(«-i)(v^-l)tfJ-(2-VS?K 9 +K 9 )j fOT g > ^/2/gK 2 . 



Therefore, we arrive at 



£ <^- 1)/f for g < y /g/gg! , (5^ 2 >2). (5.12) 
1 £ for g > ^2/gK 2 . 



These results completely reproduce those obtained so far by using various methods .EJLj'Ej 

The merit of the present approach lies in the fact t hat s patial correlations are quite easy to calculate. We assume 
that qi > q 2 without loss of generality. Noticing Eq. (3.12) leads to 



1 



\ x - y\M<lU<la) 



u<v-+n-iT(^+<iz)(y)/p{E) 



~ 1 FX? ;P (qi+q2) (E), (5.13) 

where we have kept the most relevant operator for the OPE in the second line. Some comments are in order: First, 
the OPE used her e is v alid in th e energy scale \x — y\ -C E~ Zl , where Z\ i s the dynamical exponent in Eq. ( |5.5| ). 
Next, the formula (|5.10| ) or ( |5.12| ) applies to pfa+v*) (E) in the last line of (|5.13| ). Finally, the exponent A(q 1 ,q^f is 
calculated as 

A{ qi ,q 2 )= (2 - z qi ) + (2 - z q2 ) - (2 - z qi+q . 2 ) 

2gK 2 qi q 2 for gK 2 < 2/ ( Ql + q 2 ) 2 , 

V8gK( qi + q 2 ) - gK 2 (q 2 + q 2 ) - 2 for 2/{q x + q 2 ) 2 < gK 2 < 2/q 2 , 

V8g-Kq 2 -gK 2 q 2 for 2/q 2 < gK 2 < 2/q 2 , ^> 

2 for 2/ql < gK 2 . 



The equation shows that if disorder is strong enough, the exponent A(qi,q 2 ) saturates at 2 for any q\ and q 2 . This 
saturation implies that the model is in the fully frozen phase. 

Ryu and Hatsugai0 have recently calculated the den- ( 2g for g < 1/2, 

sity correlation function of the zero energy state of the A = < 2 (^/8g — g — l) for 1/2 < g < 2, (5.15) 

Dirac fermion model, which corresponds to q± = q% = 1 [2 for 2 < g. 

and K = 1 case. Substituting these values into Eq. 

(5.14), we have It should be noted that if one compares this formula di- 

rectly with the calculation of Ryu and Hatsugai, one has 
to replace g — > g /it. The above formula seems to fit their 
calculation, and therefore we can claim that an evidence 



G 



of the strong disorder regime has actually revealed itself 
in their numerical calculation of the exponent. 



VI. SUMMARY AND DISCUSSIONS 

In this paper, we have studied the Dirac fermion model 
with the random vector field in order to calculate the 
DOS, the IPR, and their spatial correlations both in the 
strong and the weak disorder regime. We have derived 
scaling equations of one-loop order, using the OPE tech- 
niques as well as taking account of the fusion of the vertex 
operators. We have been able to reproduce the results 
known so far. Especially, we have shown that the numer- 
ical calculation of the density correlation function can be 
explained by the RG method in this paper. It turns out 
that the saturation of the exponent at 2 is manifestly due 
to the freezing transition. 

So far we have studied the random Dirac fermion in the 
continuum limit. In this case, the fermion has of course 
chiral symmetry at any realization of disorder. On the 
other hand, starting from the lattice model, e.g, the ran- 
dom hopping fermion on a square lattice_with 7r-flux per 
plaquette,cd or on a honeycomb lattice ,EJ we find, in a 
continuum limit, a Dirac fermion with two flavors due to 
the species doubling. These lattice models could be de- 
scribed by such a Dirac fermion with a random imaginary 
vector field as well as a random mass term,E3 which has 
indeed chiral symmetry. The pioneering work by GadeEij 
of the disorder systems with chiral symmetry has pre- 
dicted that the DOS diverges toward the zero energy as 
E~ x exp(— c\ ln£J|il with x = 1/2. Recently, Motrunich, 
Damle, and Huseo have pointed out that the DOS of 
such models should obey the the same scaling law but 
with x = 2/3 due to strong dispsder effects. Surpris- 
ingly, Mudry, Ryu and FurusakEJ have quite recently 
shown that such a behavior is also due to the freezing 
transition. Their method is based on the RG equation 
of the DOS, by taking account of all relevant perturba- 
tions. The KPP equation also plays a crucial role in their 
analysis. 

In this sense, the KPP equation, which was-proposed 
originally as a biological problem of the gene,E3 is a key 
equation to understand various kinds of random systems 
such_4S spin glassesjHil a directed polymers on a Gayley 
treej^j and a random hopping fermion on a lattice.Ej 

Although we cla imed in Sec. [v| that the logarithmic 
corrections in Eq. (5.2) are universal, we have neglected 
them for simplicity. It is, however, readily seen that they 
gig the logarithmic corrections, for example, to the DOS 



2-(\/sJ-i)R: 



p(E) ~ E tV»i-v«: \\ n E\ 



where 



for 
for 



gK* = 2, 

gK 2 > 2. 



(6.1) 



(6.2) 



In deriving the equation, we have included the logarith- 
mic corrections into the dynamical scaling exponents via 
the relation z q — (diut yp )q — [dim(l) — K]q. In the 
weak disorder regime, however, no logarithmic correc- 
tions should exist in the DOS, as can be seen from Eq. 
©. 

In pure systems without disorder, it is well-known 
that marginal perturbations give rise to such logarithmic 
corrections. In random systems, however, logarithmic 
corrections are allowed even in the fixed-point theories 
themselves_due to the presence of so-called logarithmic 
operators.Ej In the present case, it may be quite interest- 
ing to address the question on the origin of the logarith- 
mic corrections appearing in the KPP equation. 
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